Galaxy growth by merging in the nearby universe 
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ABSTRACT 

We measure the mass growth rate by merging for a wide range of galaxy types. 
We present the small-scale (0.014 < r < 11 h? Mpc) projected cross- correlation 
functions w(r p ) of galaxy subsamples from the spectroscopic sample of the NYU 
VAGC (5 x 10 5 galaxies of redshifts 0.03 < z < 0.15) with galaxy subsamples 
from the SDSS imaging (4 x 10 7 galaxies). We use smooth fits to de-project the 
two-dimensional functions w(r p ) to obtain smooth three-dimensional real-space 
cross-correlation functions £(r) for each of several spectroscopic subsamples with 
each of several imaging subsamples. Because close pairs are expected to merge, 
the three-space functions and dynamical evolution time estimates provide galaxy 
accretion rates. We find that the accretion onto massive blue galaxies and onto 
red galaxies is dominated by red companions, and that onto small-mass blue 
galaxies, red and blue galaxies make comparable contributions. We integrate over 
all types of companions and find that at fixed stellar mass, the total fractional 
accretion rates onto red galaxies (~ 3 h^o percent per Gyr) is greater than that 
onto blue galaxies (~ 1 h 70 percent per Gyr). These rates are almost certainly 
over-estimates because we have assumed that all close pairs merge as quickly as 
the merger time that we used. One conclusion of this work is that if the total 
growth of red galaxies from z = 1 to z = is mainly due to merging, the merger 
rates must have been higher in the past. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: 
fundamental parameters — galaxies: interactions — galaxies: general — meth- 
ods: statistical 
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INTRODUCTION 



The galaxy mergers may play an important role in the evolution of the galaxies. In 
the color-magnitude space, galaxies are separated into two distinct regions: (1 ) the 'red 
sequence': the 'early-type', red galaxies ; (2) the 'blue cloud' or 'blue sequence': the 'lat e- 



type', blue galaxies with strong ongoing star formation (IStrateva et al.ll200ll ; iBlanton et al. 



2003bl ). Some recent studies in the hig h-redshift 



1 ) universe find that the early-type 



galaxy popu 


ation is growing over time (Bell et al. ! 


2004: Willmer et al. 2006: Blanton 


2006; 


Brown et al. 


2007: Faber et al. 


2007; 


Zhu et a. 


. 201 


1), which is also ] 


bund at very hig 


i (z ~ 


2) redshift ( 


Daddi et al. 2005: 


Truiillo et al. 


2007: 


Longhetti et al. 


2007: Toft et al. 


2007; 


Conselice et al. 20071: Cimatti et al.l 


2008: van Dokkum et al.l 2008: 


Saracco et al.l 2009). In 



numerical simulations, some studies show us that major mergers of intermediate-stellar-mass 
late-type galaxies maybe play an important role in the growth of the inter mediate-stellar- 
mass early type galaxies (IBarnes fc Hernquistl Il996l ; iNaab fc Burkertl 120031). However, the 
massive early type galaxies may grow in a different way (INaab fc Burkertll2003l ) . Some recent 
studies show us that 'dry mergers' - the mergers between ear ly- type galaxies - might play an 
important role in the gro wth of massive early- type galaxies (jBell et al.l l2006bl ; Ivan Dokkum 
20051 : iMasiedi et alil2008r i. 



There are a lot of studies that estimate the merger rate among galaxies. These studies 



merger' close pairs and convert the 'pre-merger' pairs to a merger rate (e.q 


. ZeDf & Koolll989 




Carlberg et al. 


1994. 


2000: Carl 


jerg et al. 


2000; 


Patton et al 


1997 


2000: 


van Dokkum et al 


1999 


Lin et al. 


2004; 


Bell et al. 


2005 


Bell et al. 


2006a: 


. Carta 


Jtepe et al. '. 


2007: Im 


asiedi et al 


2006. 


2008:IPatton & Atfield 20081: Kitzbichler & White 


2008 


Bundv et al.ll2009: 


De Propris et al. 


2007: 


Lin et al. 


boos] 


de Ravel et al. 


2009|; 


De Propris et al. 


>010; ] 


■lobaina et al. 


2010). The 



studies in the second category count the 'post-merger' galaxies which have recently experi- 
enced at least one merger event. These 'post-merger' galaxies are chosen by some observable 
special properties caused by merging. An exam ple might be by star formation indicators 
of 'post-merger ' galaxies (IQuintero et al.l 120041) or by morphological signature s caused by 



merger events (jAbraham et al. 



20061 : iDe Propris et al 



2007 



1996; 



Conselice et al 



Lotz et al 



2003; 



20081 : IConselice et alihool? ) 



van Dokkum 2005: Lotz et al. 



Our work builds on the earlier works of iMasjedi et al.l (120061 . 120081 ). which have found 
previously that luminous red galaxies (LRGs) are gro wing on average by l ess than 2 h-? percent per Gyr 
from merger activity at redshifts 0.16 < z < 0.30 (IMasjedi et al.l 120081 ) . In this paper, we 
consi der both red and blu e galaxies. We use the previous technique for measuring the close 
pairs (IMasjedi et al.ll2008l ) on NYU VAGC Spectroscopic Sample and SDSS Imaging Sample, 
and extend this type of analysis beyond luminous red galaxies (LRGs) to a wide range of 
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galaxies in both stellar mass and color. 



The primary uncertainty in turning a deprojected three-dimensional cross- correlation 
function at small scales into a merger rate is in estimating the mean time for two galaxies to 
merge as a function of stellar mass and separation. There are different estimates of merger 
timescale: free-fall time, orbital t i me, and dynamical friction time (e.g.. IBinnev fc Tremaine 



2009 



1987; Bo vlan-Ko: 



Lotz et al 



chin fc Mall2007l: IConrov et al1l2007l : iKitzbichler fc WhitelboOsIlBundv et al 



20101 ) . In this paper , we will use an app r oxima tion to the two kinds of 



merge r times under the assumption of IKitzbichler fc White! (120081 ) and IBinnev fc Tremaine 
(119871 ) as our standard estimate. Both of these times are likely to be an under-estimate of 
the mean merger time, because some close pairs will not merge at all. Any under-estimate 
of the merger time leads to an over-estimate of the growth rate. 

Throughout this paper, all magnitudes are AB, all apparent magnitudes are model 
Mag, all masses are stellar masses (in units of /if 2 M Q ), all velocities are in u nits of km s~ 



all rad i i of galaxies are rqn which c ontain 90 percent of the Petrosian flux (IBlanton et al. 



2003bl ; IBlanton fc Moustakad 120091 ). and all volumes and distances are comoving, calcu- 
lated for a cosmological world model with (Q m , Q\) = (0.3, 0.7) and Hubble constant H = 
70h 70 kms^Mpc- 1 . 



2. DATA 



We use the NYU Value-Added Galaxy Catalog (VAGC) V7.2 data (IBlanton et al. 



20051 ). which is built from SDSS data, as our source of spectroscopic data. We use all 



vey of about 10 4 square degrees ( 


Fukurita et al. 


1996: 


Gunn et al. 


1998 




Gunn et al. 


2006; 


York et al. 


1999: 


York k SI 


DSS Collaboration 


2001: 


Stouehton et al 


2002 




Abazajian et al. 


2003; 


Abazajian et al. 


2004 


)■ 







2.1. Spectroscopic Subsample 



Our spectroscopic sample is drawn from NYU VAGC V7.2 data. We removed galaxies 
with apparent magnitude s r < 14 mag, b ecause SDSS photometric catalog missed many 
luminous galaxies nearby ( IZhu et al.l 120101 ). Our spectroscopic sample c ontains about 8.6 x 
10 5 main galaxies, of which about 5 x 10 5 SDSS Main Samples galaxies ( jStrauss et al.l 120021 ) 
are in the redshift range 0.03 < z < 0.15 (see Figured]) with apparent magnitudes r < 17.77. 
We cut the redshift at z < 0.15, because we want to avoid all quasars, and we cut the 
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redshift at z > 0.03, bec ause we want to avoid all stars. We use the code sdss_kcorrect 
(IBlanton fc Roweid 120071 ) to calculate the if-corrected absolute magnitude at z = 0.1 and 
stellar mass for each. We cut this sample into 10 subsamples by stellar mass and make sure 
that in each subsamples th e number of galax ies is the same. After that, we cut them into 
red and blue using the cut (IHogg et al.l 120041 ): 



o.i 



9 



o.i 



0.0625 x log(M) + 0.15 



where M is the stellar mass. These cut the spectroscopic sample into 20 spectroscopic sub- 
samples; we call these subsamples D s (see Figure [2]). For the red spectroscopic subsamples, 
we name the smallest stellar- mass subsample 'redO', then name the second smallest stellar- 
mass subsample 'redl', and so on, so the largest stellar-mass subsample is named 'red9', 
and the same for the blue spectroscopic subsamples. Figures [3] and H] show the number den- 
sity of galaxies of the 18 spectroscopic subsamples as a function of redshift. The redO and 
blueO subsamples are low in number density, and highly affected by survey selection effects, 
so we do not use them further. We can see that blue8, blue9 and red9 appear to rise in 
number density with redshift, this is because we removed galaxies with apparent magnitudes 
r < 14 mag. 

The SDSS suffers from the 'fiber collision': the angular separation between any two 
spectroscopic targets must be larger than 55 arcsec. There are about ~ 9 percent of target 
galaxies which do not have redshi f ts bec ause of this fiber collision. Using the counting-close- 
pairs technique of iMasjedi et al.l ( 120081 ). the fiber collision will not affect our pair counts 
directly under the approximation that the unmeasured galaxies and the measured galaxies 
are similar in cross-correlation with fainter galaxies. 



Following the technique of lMasjedi et al.l (120081 ) . we estimate the weight pj that accounts 
the spectroscopic incompleteness from fiber collisions effects, and the weight fj that accounts 
the spectroscopic incompleteness from all other selection effects in SDSS. For each galaxy j 
in the spectroscopic subsample, we calculate pj by using a two-dimensional 'FOF' (friends- 
of-friends) grouping algorithm on our main galaxies targets with a 55-arcsec linking length: 



Pj 



JV total 



(2) 



Where N 3 total is the total number of main galaxies targets in group j, N-> is the number of 
main galaxies with redshift measurement in group j, and the group which contains galaxy 
j, is called group j. All pj > 1 because N^ otal > iV|. 

For each spectroscopic subsample D s , we create a random spectroscopic subsample M. s 
which is 10 times large as D s . For each galaxy in D s , we created 10 galaxies with the exactly 
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same redshift but angular position taken from the two-dimensional random sample. Thus, 
we create a large random spectroscopic subsample R s which match the redshift distribution 
of D s , and R s is isotropic within the SDSS survey region. 

The SDSS survey region is separated into small unique region 'sector'. For each random 
point j in R s : 



fj = — (3) 



1 

Where Fj is the fraction of main galaxy targets for which a classification was obtained in 
the object's sector (in our NYU-VAGC spectroscopic sample, the average Fj is F ~ 0.91). 



Our correction of fiber collisions is: we weight the target j in spectroscopic subsample D s 
as pj and weight the target j in random spectroscopic subsample R s as fj . From the previous 



work (jMasjedi et al.ll2006l ). we know that this correction will improve in the spectroscopic 
incompleteness due to fiber collisions at very small sep arations (w(r p ) < l OOhjn kpc). In 



Section |4"~T| we will compare our result with the result of IZehavi et al.l (120101 ). and will show 
that our result with this correction of fiber collisions fits better than the result without this 
correction of fiber collisions. 

For each spectroscopic sample 's' (D s or R s ), in which there are N s spectroscopic galaxies, 
we divide it into the 50 bins by lines of constant dec, so that there are iV s /50 spectroscopic 
galaxies in each bin. Then we resample them into 50 leave-one-out resampling samples, 
so that there are 49/50 x N s spectroscopic galaxies in each sample. We call them the '50 
resampling samples', with which we can calculate our jackknife resampling covariance matrix. 



2.2. Imaging Subsamples 

For our imaging data, we use a sample drawn from the full SDSS imaging catalog in 
which there are about 4 x 10 7 galaxies. We include from the SDSS imaging sample only 
galaxies with apparent magnitude 14 < r < 21.5 mag, and apparent color —0.5 < [g — r] < 
2 mag, see Figures [5] and [61 We removed galaxies with apparent magnitudes r < 14 mag for 
the same reason with our spectroscopic sample. We also removed galaxies with apparent 
magnitudes r > 21.5 mag, because these galaxies are not well observed and their observed 
number density is much lower than their real number density. Please note that this cut will 
affect our minor merger near mass ratio 1 : 100. Similar to the spectroscopic sample, we 
create a random imaging sample Rj as large as possible. The angular positions of galaxies in 
Rj are taken from the two-dimensional random sample. So R« is isotropic within the SDSS 
survey region. 
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2.3. Grid Method ^-corrections for galaxies from the imaging sample 

We also need to compute the stellar mass of galaxies in the imaging subsamples in order 
to determine the mass ratio, so we want to i^-correct galaxies in the imaging subsamples. 
However, we cannot f^-correct individual galaxies in the imaging subsamples once and for 
all, because we do not have spectroscopic redshifts for them. Each time we consider a pair 
of galaxies, one from the spectroscopic subsample and one from the imaging subsample, we 
assign the spectroscopic redshift to the galaxy from the imaging sample. This allows us 
to calculate for each galaxy from the imaging sample in each spectroscopic-imaging pair a 
temporary .fT-corrected stellar mass and [ 0A g — 01 r] color for the purposes of that pair. We 
discard these values and compute new ones when the galaxy from the imaging sample is used 
in another pair with another galaxy from the spectroscopic sample. 

To save time, we take galaxies from the NYU VAGC spectroscopic sample as represen- 
tative of all galaxy types, and apply the code sdss_kcorrect on the galaxies from a grid 
named 'B' of observed r-band magnitude (0.5 mag per bin), [g — r] color (0.1 mag per bin), 
and redshifts between 0.03 and 0.15 (0.0002 per bin). We saved the mean i^-corrected stellar 
mass Mb and [ ai <?,B — 01 rs] color in a grid of observed r-band magnitude, [g — r] color, and 
redshift, also save the mean redshift zb, mean tb and mean [g B — r#] color, Thereafter we 
estimated the i^-corrected stellar mass and [ 0A g — 01 r] color for a galaxy G in grid B, and 
the [go — r G ] color of galaxy G is between [g B — r B ] color of grid B and [gc — r c ] color of 
grid C which is next to grid B (that means zb and z G are the same, and r# and r G are the 
same): 

, M , fM v (d G L )\ r G -r B (log M c - log Mb) x {[g G - r G ] - [g B - r B \) 

log M G = log {M B x t-tbt^ ^7- + f j f j 

[dlY 2.5 [g c - rc] - [g B - r B \ 

(4) 

[ 0A 9G ~ 01 r G ] = [° A 9b - 0A r B ] + [g G - r G ) - [g B - r B ] (5) 

where is the luminosity distance calculated from redshift, and galaxy G is in grid B, so 
the difference between zb and z G is small, similarly, the difference is small between d'f and 
c?f, Tb and r G , [gs — re] and [g G — r G ]. because of these small differences, the equations 
(jl]) and (ED can be used. This speeds up the correct ion procedure immensely and only 
introduces a 12-percent one-sigma error for each galaxy and there is little bias (see Figure [7]), 
so only introduces percent-level errors in the results. We call this 'Grid Method' hereafter. 

For grids with observed r-band magnitude r > 17.77 mag, we cannot get the mean stellar 
mass and [ 01 g — 01 r] color directly, because there are no galaxies from the spectroscopic 
sample at observed r-band magnitude r > 17.77 mag. In order to estimate the color and 
stellar mass of the galaxies in grid A with mean observed r-band magnitude ta > 17.77 mag, 
we find the grid point B that has the nearest mean observed r-band magnitude r# (r# > 



-7- 



17.77 mag) and has the same observed [g — r] color and the same redshift, also there are 
at least 10 galaxies from the spectroscopic sample in grid B. Then we estimate Mq and 
[ 0A 9g — °' lr c] color of galaxy G in grid A using the equation (jl]) and (J5J). If we cannot find 
a grid B satisfying the conditions, we will leave all the galaxies in grid A empty 

Figure [8] shows the difference between the color estimate using the Grid Method above 
and the color calculated using the code sdss_kcorrect. The one-sigma error for each galaxy 
is about 0.04 mag and there is little bias (see Figure [8]). 



2.4. Velocities for merger rate estimates 

In this Section, we will estimate the average orbital velocity which will be used to 
estimate the merger time in Section 13.41 Please note that this approximation will induce a 
large error in the estimate of merger time, because we assume that all close pairs merge under 
the following orbital velocity which is not true for close pairs in high velocity-dispersion. 

The average orbital velocity for a galaxy from the imaging sample around a more massive 
red galaxy from the spectroscopic subsample s with average velocity dispersion a v is very 
roughly 1.5 t imes the velocity d ispersion, here we have included the factor of 1.5 to be 



conservative (jMasjedi et al.ll2008l ). We estimate the a v with Faber- Jackson relation: 

log a v = ai + h log M r s ed (6) 

where M r s ed is the mean stellar mass of the red galaxies from spectroscopic subsample s. We 
performed a linear fit to the data to obtain a\ = —1.588 and b\ = 0.354 (Figure [9]). Please 
note that the method above will induce a small enough error (< 1 percent) into our final 
result. 

We estimate the average orbital velocity V c for a galaxy from the imaging sample around 
a more massive blue galaxy from the spectroscopic subsample s with the Tully-Fisher rela- 
tion: 

log V c = a 2 + b 2 log (Li) (7) 

where (Lj) is the mean I-band luminosity of the blue galaxies from spectroscopic subsample 
s, calculated from Lj = M/rj, where M is the stellar mass of the galaxy and rj is the I-band 
mass-to-light ratio of the galaxy calculat ed from the code sdss _kcorrect. We used for this 



relationship a 2 = —0.835 and b 2 = 0.291 (ICourteau et al.l 120071 ) 



Table Q] provides this information for all 20 spectroscopic subsamples. 
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Method 



For each spectroscopic subsample D s , we cut the imaging sample into 16 subsamples 
by stellar mass: 10~i < Mf '/M 8 < 10" V, where Mf is the stellar mass of a galaxy from 
the imaging sample, calculated by the Grid Method using the redshift of the spectroscopic 
galaxy, and M s is the mean stellar mass of the galaxies from spectroscopic subsample s, and 
j is an integer 1 < j < 16, that means a mass ratio of 1 : 1 to 1 : 100 is covered. After 
that, we cut them into red and blue using equation (QQ) by [ 0A g — 01 r] color calculated by 
Grid Method. For each spectroscopic subsample D s , these cut the imaging sample into 32 
imaging subsamples Oj. 

In this section, we will show our method to estimate the merger rate between galax- 
ies in spectroscopic subsample D s and galaxies in imaging subsample D^: (1) we estimate 
the projected two-dimensional cross-correlation function w 



siy Pj 



as a function of tangential 
projected separation r p ; (2) we de-project the smooth fit for the cross-correlation function 
w S i(r p ) to obtain the three-dimensional real-space cross-correlation function £ a »(r) as a func- 
tion of real-space sepa ration r; (3) we estimate t he merger rate using £ S j(r ) and our two 
kinds of merger times (IBinney fc Tremaind Il987t iKitzbichler &: White! 120081 ) . In addition, 
we will also discuss our method of photometry correction. 



3.1. Projected cross-correlation function 



To estimate the w si {r p ) between spectroscopic subsample D, s and imag ing subsample D, 



we can integrate £ S i(r) along the line of sight {e.g., iDavis fc Peebles! Il983h : 

poo 

w sl {r p ) =2 dy£ si [(rj + y 
Jo 

This integral is dominated by scales y < r p . 



, 2 \ 1/21 



Using the previous approach (IMasjedi et al.l 120061 . 120081 ) . we estimate not w 



J siy pj 



but 



PiW S i{r p ), where pi is the average comoving three-dimensional stellar mass density of the 
imaging subsample v. 

D s Di R s Di 

Piw si {r p ) = — — - — — (9) 
U s Ki rC s Ki 

Where D s and Di represent the spectroscopic and imaging data subsamples, and R s and 
Ri repre sent the spectroscopic an d imaging random subsamples. Similar to the previous 



method flMasiedi et al.ll2006 



20081 ) . equation Q measures the mass- weighted abundance of 



pairs (j H^) and subtracts th e mean background level This method has been well 



tested in Masiedi et al.l (120061 ) 
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In detail, the factors are defined as follows: 



Pi 



D s Di 



where the top sum counts the weighted pairs j of galaxies from D s and ID 
tangential projected distance r p , and the bottom sum is over galaxies j from 
D s Di is dimensionless. 



(10) 



separated by 
This factor 



D s Ri 



an 

dA 



J J 



dM 
dfi 



(11) 



where the top sum counts the weighted pairs j of galaxies from D s and Mj separated by 

) is the inverse square o 

5 S , and 4^ is calculated by 



tangential projected distance r p . In the bottom, (^) is the inverse square of the transverse 

is calculated by - dN 



comoving distance (lHoggill999l ) to galaxy j from 

dN 
dfl 



where 4tt is the two-dimensional number density of Rj per solid angle, and 



121 



where M{ is the mean stellar mass of galaxies from the imaging subsample i, and j is an 
integer 1 < j < 16. Then (j^j x ^ represents the average stellar mass of galaxies in Rj 

per unit comoving area around each galaxy from D s . This factor -D s i?i has dimensions of 
comoving area divided by stellar mass. 



£ u 



£a 



(13) 



this is similar to equation ( FlUj) . but M s represents now the random catalog mentioned in 
Section 12.11 This factor R s Di is dimensionless. 



E u 



R s Ri 



jel 



dA 



dM 

~dn 



(14) 



this is similar to equation ( II ip . but R a represents now the random catalog mentioned in 
Section 12.11 This factor R s Ri has dimensions of comoving area divided by stellar mass. 
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For some experiments, we need to estimate niW S i(r p ), where is the average co moving 
three-dimensional number density of the imaging subsample i. We estimate this by the 
following estimation: 

/ \ D S D. L R s Di 
rii w si {r p ) = - (15) 

[L> s Ki\ N [K s Ki\ N 

For D s Di and R s Di, it is as the same as equation fflQ|) and equation ffT3]) . For [D s Ri]N 
and [R s Ri\]y: 



[^^]iV = r.nl ( 16 ) 



dfi 
dl 



J j 



dN 

dn 



this is similar to equation f TTTj) . but ^ is the two-dimensional number density of the random 
imaging catalog per solid angle. This factor D s Ri has dimensions of comoving area. 

E h 

\RM» = — ^ JT, (17) 



cL4 



dN 

dn 



this is similar to equation f fl4|) . This factor -R s -Rj has dimensions of comoving area. 

In range of our interest 0.0149 < r p < 11.9 h^ Mpc, we bin the spectroscopic-imaging 
pairs counting by the comoving projected separation r p of the pair where r p = = 0.0149 x 
10 fc / 5 /ifo 1 Mpc and k is an integer < k < 14. We have already discussed how to bin the 
spectroscopic sample D s and the imaging sample D» in Section [2j We have combined the 16 
stellar mass bins into 4 to simplify the figures. Figures [TD] through [T5] show the results of 
our measurements of PiW si (r p ). 

The uncertainties on the results shown in these figures are estimated using jackknife 
resampling covariance matrix with 50 resampling samples (see Section [2TT|) . please note that 
all the error bars in our graphs only come from the jackknife resampling covariance matrix 
(there are some other errors like the error of the color and stellar mass estimated by the Grid 
Method, and so on). On hundreds of kiloparsec scales, the error bars for each subsample 
are smallest. On smaller scales, the error bars become larger because of the 'shot noise': 
the smaller the separations, the fewer the pair counts. On larger scales (> 1 Mpc), the 
error bars become larger because there are more and more interlopers on larger scales which 
means that the background subtraction is more noisy. 

Figures [TUl through [131 show that niW S i(r p ) is a complex function of r p . However, on 
very small scales - tens of kiloparsec scales, we assume that rij w S i(r p ) scales (something) like 
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r 1 (IMasjedi et al.ll2006l . 120081 ). We fit each set of w S i(rk) data with the smooth model: 



w(r p ) = w 



w 


1 + 


7 


V 











-1 



by minimizing 



(Piw(r k ) - pjW si (r k )f 



(19) 



We choose r c = 12.5 rgo for Figure [10] to Figure [13] (where rgo is the median radii of the 
galaxies in the corresponding spectroscopic subsamples) ; then find the 7 and pi w that 
minimizes xli- Figures IT0l through [TBI show us these fits. 



3.2. Photometry Correction 



One important issue with all clustering measurements on small scales is possible pho- 
tometric biases when measuring cl ose pairs. This issue can directly lead to biased flux 
measurements (IMasjedi et al.l 120061 ) and biased color measurements for galaxies, and will 
indirectly affect s tellar masses, k-cor rections, etc. This can be due to poor photometry in 
crowded systems (IPatton et al.ll201ll ). 



We b uild our method to cor rect photometric biases upon the photometry test of lMasjedi et al. 



(120061 ). In lMasjedi et al.l ( 120061 ). they created fake images of pairs of identical galaxies with 
separations ranging from 2 to 35 arcsec. These galaxies represent passively evolving LRG 
galaxies observed at a redshift of z = 0.3 with de Vaucouleurs profiles (n = 4 Sersic profiles). 
Then they placed one such galaxy pair onto RUN 2662 (which has a typical SDSS seeing 
of about 1 arcsec) of SDSS imaging. After inputting the known info into the mock galaxy 
images, they processed these images as raw SDSS images using the standard SDSS pipeline, 
PHOTO, to determine the effect of proximity of galaxies on their measured properties (see 
plot 1 in Figure [T41) . At separations larger than 20 arcsec, the Petrosian flux measures 79.5 
percent of the input Sersic flux, which is calculated by three sigma out-layer rejected aver- 
age. In other words, the Petrosian flux only measures about 80 percent of a galaxy's light. 
We are interested in intermediate separations, (5 < s < 20 arcsec), in which the fraction 
of the recovered flux to input flux increases to 83 percent. This increase is likely due to a 
double counting of the low level diffuse emission from the two galaxies which is being poorly 
deblended between the two objects. 

For pairs of main galaxies, we study two different cases: one for galaxy pairs consisting 
of two identical galaxies and another with galaxies of different stellar mass. For Case 1, we 
consider a pair of identical main galaxies with radii of r^ am (at redshift 0.03 < z main < 0.15) 
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than that of the LRGs r^ RG at redshift z = 0.3 (see plot 1 and 2 in Figure [T6|) . We take 
the following two approximations: (1 ) If we consider a pair of identical main galaxies which 
are at the same redshift z = 0.3, then the only difference between this pair of main galaxies 
and the pair of LRGs (see plot 1 in Figure [T4"j) is that the angular radii of the main galaxies 
is smaller than that of LRGs by a factor of r™ in /rg RG . So we compress the result of plot 1 
in Figure [TH by a factor of r^ am /rg (see plot 2 in Figure IT41 as an example) as the effect 
of proximity of a pair of identical main galaxies. (2) If we consider a pair of LRGs which 
are at a different redshift z = z ma i n , then the only difference between this pair of LRGs 
and the pair of the LRGs at redshift z = 0.3 (see plot 1 in Figure [T4"j) is that the angular 
radii of this pair of LRGs is larger than that of the LRGs at redshift z = 0.3 by a factor 
of D mam I D LRG , where D main is the comoving distance from z = z ma i n to us, and D LRG is 
the comoving distance from z = 0.3 to us. So we stretch the result of plot 1 in Figure HH 
by a factor of D mam / D LRG (see plot 1 in Figure US as an example, and D mam at redshift 
z = 0.144246 is nearly half of _D ii?G ) as the effect of proximity of a pair of LRGs at redshift 
z = z ma i n . Combining the above two approximations, the final effect of proximity of a pair 
of identical main galaxies at a different redshift z = z ma in 7^ 0.3, will be stretched by a factor 
of r™Q %n /r$ RG x D mam j jjlrg ^ gee 2 Figure [TBI as an example). Using this method we 
correct the flux measurement of our sample on small scales of major merger between main 
galaxies and we assume that the correction of stellar mass is equal to that of flux. For Case 
2, it is based on the above Case 1 but involving the radii rgo of the pair of galaxies (ta and 
Tb) in our estimation. In plot 2 of Figure [T6| we take the third approximation: (3) The flux 
density from the left galaxy onto the right galaxy is a constant Da (this approximation lead 
to a small error comparing to the flux density from the galaxy A onto the galaxy B in plot 
3 of Figure [IB]) . We mark the percent increase of flux from the other galaxy as P/iux, that 
of stellar mass as P maS s, radius of galaxy A as ta total flux as Fa and stellar mass as Ma- 
Then we get the following result: 

Da x 7r r\ 

Mnass * flux ^ 
?A 

In plot 3 of Figure [161 we assume that the flux density from the galaxy A onto the galaxy B 
is equal to Da and the flux density from galaxy B onto galaxy A is equal to Dg. We mark 
the percent increase of flux of galaxy B from the galaxy A as Pfi ux , that of stellar mass as 
Pm ass , radius of galaxy B as r# total flux as Fb, stellar mass as Mb, and the percent increase 
of flux of galaxy A from the galaxy B as Pfi ux , that of stellar mass as Pm ass . Then we get 
the following results: 

Db = DaX- = D a x— (21) 
r>A r>A Db x 7rr^ M B 

"mass = *flux = ~p = *mass x l^Zj 
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B B Da x 7r r| M A r% 

"mass = flux = 7? = *rnass x ~T7 X ~ \ ^/ 

Using this result we correct the flux and stellar mass measurements of our sample on small 
scales of minor merger between main galaxies. 

Please note that our photometry correction is overestimated because of the above ap- 
proximations (1) and (2). In approximations (1) and (2), we assume that the absolute 
angular scale does not matter; and the only change comes from the ratio of the absolute 
angular radii of the pair of galaxies to the absolute angular separation of the pair of galaxies. 
However, we know that absolute angular scale does matter: at the same ratio the larger the 
absolute angular is, the easier the deblending will be. So approximations (1) and (2) will 
contribute a few percent error in our final result of photometry correction. 

After the above photometry correction, we reset our spectroscopic and imaging subsam- 
ples to recalculate the PiW S i(r p ) using the method in Section [37T1 (see the first and second 
data points in Figures [TU1 through [T31) . Please note that we also apply our photometry cor- 
rection on the galaxies with no nearby companions, however these galaxies have zero weight 
in our pair-counting and will not affect our result, because there are no companions near 
these galaxies during counting pairs. In order to show our method of photometry correction 
is robust, we double our photometry correction and find that all the percentage difference 
between the result of our photometry correction and double our photometry correction is 
below 26 percent for one data point. Then, this data point with 26 percent change will only 
contribute a few percent error in our final result of the total fractional accretion rate after 
our fitting curve (see equation (1TB]) and (|T9|) ). So, if we assume that the percent error of our 
photometry correction in flux is 100%, the final effect onto the total fractional accretion rate 
is at most a few percent. 

We use the above method to correct the photometric bia ses, and we find that o ur correc- 



tion due to photometric biases is much smaller than that of IMasjedi et al.l ( 120061 ). because: 



(1) Main galaxies have smaller radii than LRGs, so it is easier to deblend a pair of main 



galaxi es than a pair of LRGs. (2) Photometry correction of auto-correlation of lMasjedi et al. 



( 120061 ) is larger than that of our cross-correlation, this is because of the difference of the stellar 
mass cut of spectroscopic/imaging sample between us: for auto-correlation, the stellar mass 
cut of spectroscopic/imaging sample will be M s > M l t ™^ hold and Mj > M^^ hold , so after 
photometry correction, the only effect is that some spectroscopic (and imaging) galaxies near 
^threshold will be cut off from the spectroscopic (and imaging) sample which will decrease the 
pair-counting. However this is not the only effect on cross-correlation. For cross-correlation, 
the stellar mass cut of spectroscopic and imaging sample will be M l t ^!^ hold < M s < M^^ hold 
and M l ° wer < Mj < M^ pper , so after photometry correction, besides the above effect there 
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is another effect that some spectroscopic and imaging galaxies which are a little bit above 
^threshold or M^ pper will be counted into the spectroscopic/imaging sample from outside. This 
effect will increase the pair-counting. Combining the two effects above for cross-correlation, 
the final photometry correction for cross- correlation will be smaller than that for auto- 
correlation. 



3.3. Three-dimensional statistics 

The smooth fit w{r p ) to each projected correlation functions w si {r p ) can be deprojected 
to get an estimate of the three-dimensional space correlation function £(r) by 

M(r) = ~ rdr p d[pi 2 irp)] (rl-rr l/2 (24) 



7T 



e.g. 



Davis Peebles! Il983l ). where pi is a constant. 



The mean total stellar mass M * of galaxies from a specific imaging subsample i within 
a given small three-dimensional separation r c / ose around each galaxy from spectroscopic sub- 
sample s is: 

M* = p t [ dV t [1 + £ si (r)] = 4vr Pi f ^Vdr [1 + £ si (r)] (25) 



o 



At small scales, £ S j(r) 3> 1, so: 



M*^4vr / r 2 dr[ Pi ^ si {r)} (26) 



o 



From [p;£si( r )] we can see that we do not need to measure pi and w S i{r p ) separately. 
At very small scale (r c / ose <C r c ): 

M* tti Pi w r c r close (27) 

Similarly, the average number N* of galaxies from a specific imaging subsample i within 
a given small three-dimensional separation r c / ose per galaxy from spectroscopic subsample s 
at very small scale {r c i ose <C r c ) is: 



N* » 4 m w r c r dose 



(28) 
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3.4. Merger rate 



We can estimate the merger rate Tj of galaxies from sample i into galaxies from sample 
s per spectroscopic galaxy per unit time by: 

NT 



f , 



(29) 

u merge,i 

The mean fractional stellar-mass accretion rate of galaxies from spectroscopic subsample s 
from merging with galaxies from imaging subsample i per unit time is: 

M* 



"dlnilC 


1 


~dM s - 




,• ~ M~s 





tmerge,i 



(30) 



In principle, all merger rate estimates depend on the radius r c i ose inside of which we 
have counted close pairs. However in this work, we are interested in the instant merger rate 
estimates, which means that we are interested in the range of r c \ ose <C r c , and this range 
can be reached using our fit lines. Another reason why we use our fit lines instead of using 
the data points at very small scales, is that the error bars due to the shot noise for each 
subsample are very large. From these fit lines in Figures [10] through [13] we know, over the 
range of interest (r c / ose <C r c ) u> S j(r p ) scales like r" 1 , £(r) scales like r~ 2 , N* and M* scale 
(something) like r c i ose . Similarly, both of the time-scales {txw,i and tBT,i) scale like r c / ose . 
For this reason, at r c i ose <C r c the above merger and accretion rates do not depend strongly 
on r dose . 



In this work, we u se these two merger ti me es timates: tsT,i from iBinney &: Tremaine 



( 119871 ) and txw,i from iKitzbichler &: White! ( 120081 ). Both of them depend on the orbital 
merger time t or ut- 



t 



red 
orbit 



2 7T r dos , 
1.5 a v 



(31) 



Where 1.5 a v is the average orbital velocity for a galaxy from the imaging sample orbiting 
a more massive red galaxy from the spectroscopic sample with velocity dispersion a v , see 
equation ([6]). Similarly, the orbital merger time for a galaxy from the imaging sample merged 
into a more massive blue galaxy from the spectroscopic subsample is 

2 7T Tclose 



t 



blue 
orbit 



V, 



(32) 



Where V c is the average orbital velocity for a galaxy from the imaging sample around a more 
massive blue galaxy from the spectroscopic sample, see equation (J6j). 



For the assumption of txw,i from IKitzbichler &: White! (120081 ). the approximation be- 
comes 

tKW.i t or f)it 



0.3 



(33) 
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where we assume M s > Mi. The solid lines in Figures [T71 and [181 as a function of th e mass 
ratio Mi/M s show the merger rate under assumption of IKitzbichler &: White! (120081 ). The 
total fractional accretion rate is the area under each curve. 



For the assumption of £bt,i from iBinney fc Tremaind (119871 ). tsT,i is long er than the 
rbita l merger time t nr kit by a factor roughly equal to the ratio of the stellar masses (IBinney &: Tremaine 
19871 : iMasiedi et alfeooa ): 



t 



BTa 



t 



orbit 



Ms 
Mi 



(34) 



where we assumed M s > Mi. The dashed lines in Figures [T7] and [18] as a function of the 
mass ratio Mi/M s show the merger rate under assumption of all possible mergers takin g 
place within one dynamical friction time (IBinney fc Tremaind 119871 ; iMasjedi et al.l 120081 ). 
and the total fractional accretion rate is the area under the curves. 

We integrate the mass accretion rate from mergers over all imaging subsamples and find 
the total fractional accretion rate (hjo Gyr -1 ) of all the main galaxies. This is shown with 
solid lines (for t KWji ) and dashed lines (for tsT,i) in Figure [T^l Please note that the errors 
from merger time are shown in this figure. 

Table [2] shows the total fractional accretion rates of all the 20 spectroscopic subsamples 
for both the assumption of t merge ^ = txw,i and the assumption of t merge ^ = tBT,i- 



4. Comparison to Previous Work 

4.1. Comparison to previous clustering results 

We estimate not w p {r p ) but niW p {r p ). In order to estimate w p {r p ) and compare with 
previous results, we need to estimate n,: 

n iWp {r p ) 

w P {r P ) = (35) 

Tli 

It is difficult within the SDSS data to precisely measure the real-space number densities 
for the imaging subsamples with stellar masses Mi < 6 x lO 9 /if o 2 M , because there is only 
good spectroscopic information about bright members of the imaging sample. However, for 
galaxies with stellar masses Mi > 6 x 10 9 /if 2 M , the real-space number density rij for the 
imaging subsample i is measurable: 

m = f (36) 

where n| is the average real-space number density for the corresponding spectroscopic sub- 
sample s within its volume limit, and F is the mean fraction of Main targets for which 
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a classification was obtained in the object's sector, for NYU-VAGC spectroscopic sample 
F « 0.91. We assume that rii is non-evolving over the redshift range of interest. We find 
the lower and upper redshift limits of the volume-limit for the corresponding spectroscopic 
subsample s in order to calculate n|. 



Galaxy clustering has been measured at intermediate and small scales (IZehavi et al. 



20051 . l20ld : Masiedi et al.ll2006l : IChenll2009l: iLi fe Whitell2010l : White et alJl201lh . Our results 



are consistent with the results of IZehavi et al.1 (120101 ) (see Figure [2D]) • In order to generate 
our mass-threshold samples which are nearly the same as their luminosity-threshold samples, 
we calculated M*, the mean mass of the galaxies nearby their luminosity-threshold M r as 
our mass-threshold, and we cut M s > M* and Mi > M* to generate the corresponding 
mass-threshold samples. This turns our cross-correlation into an auto-correlation. In order 
to calculate nf for the two subsamples with M r < —18.0 and M r < —18.5 in Figure |2"U| we 
use the peak real-space number densities instead of the average real-space number densities, 
because there are no obvious volume-limit in these two subsamples. 



Similarly, we can measure pi and w S i(r p ) separately instead of piW 



In in order to show that our result successfully correct ed the fiber collision s, we compare 
our result with the extension of the best-fit power law from lZehavi et al.l (120101 ). see Figure 1211 
and [221 Th e exten sion dashed line is from the power fit of the first six data points of 



Zehavi et all (120101 ) in the range 0.25 Mpc < r p < 2.5 Mpc. We cut at r p » 
2.5 Mpc because there is a sharp break at r p « 2.5 ft, 70 Mpc which will be discussed at 
the end of Section [5j On the other side, this extension is very robust. In Figure [221 the 
difference is very small amo ng the three extensi on dashed lines using the first five, first six 
and first seven data points of IZehavi et al.l (j2010l ). Our result (the triangle data points) with 



correction of fiber collisions fits better than the result assuming pj = 1 and fj = 1 (the 
diamond data points). We also show the data point before photometry correction Figure 



4.2. Comparison to previous merger rate results 

Please note that the merger time error will be shown in the error bars in the figures 
from now on. 

Our results are consistent with recent measurements of the merger rates based on 



counts of close pai r s dMasiedi et al.l 



Bundy et al.l l2009i ; iDe Propris et al 



2008 



2010 



Patton fc Atfield 2008; Kitzbichler fc White 2008 



Robaina et al.l 120101 ). The low dry merger rate 



an upper limit of 1.8 h 70 percent per Gyr for massive red galaxies (red9) under assump- 



tion of t. 



merge 



tsTi here is in good agreement with a number of other estimates: at 
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z < 0.36. iMasjedi et all ( 120081 ) obtained an upper limit of 1.2 h^o percent per Gyr (converted 
from 1.7 h percent per Gyr) for the dry merger rate of SDSS LRGs with Mj < —22.75; 
at 0.45 < z < 0.65, iDe Propris et all (120101 ) determined a 5a upper limit to the dry 
merger rate of 0.7 percent per Gyr (converted from 1.0 h percent per Gyr) for galaxies 
with -23 < M(r) k+e)Z=0 . 2 + 5\ogh < -21.5 in the 2dF-SDSS LRG and QSO (2SLAQ) 
redshift survey. 



Robaina et al.l (120101 ) found that the fraction of galaxies (M > 5 x 10 10 /i 70 2 M ) in pairs 



separated between 15 and 30 kpc in 3D space is /^pf^ 70 kp ° = 0.01 at z = 0.1, which is 
calculate by (1.0 — 0.3) x F(z), because they find that 30 ~ 40 percent of galaxies in close 
pairs have r < 15 hjQ kpc separations. They also expect most of the mergers to be majors; 
i.e., with mass ratio between 1:1 and 1:4. Our result of f^ D pair 7 ° kP ° * s 0-02 at 0.03 < z < 0.15 



with m ass ratio between 1:1 and 1:4, which is consistent with the result of iRobaina et al. 

(hoioh . 



With t 



merge 



txw,i determined, we compute the volumetric merger rate (the number 
of mergers per unit time per unit comoving volume) as a function of the stellar mass of the 
primary or host galaxy. We call this the merger rate mass function (merger rate MF) and 
denote it using the variable, Figure |2"31 shows the co mparison of our majo r merger rate MF 
at 0.03 < z < 0.15 with the major merger rate MF of iBundv et all (120091 ) at 0.4 < z < 0.7 
at mass ratio m/M > 0.25. We can see that both of Figures [231 and [241 show that our results 
are consistent with these previous results. 

Our results are also consistent w ith recent merger rates predicted in the ories of galaxy 
formation in a cosmological context ([Mailer et al.ll2006l : IStewart et all 120091 ) . We estimate 
our merger rate R mg at a certain mass ratio m/M by integrating Tj in equation (1291) . We 
compare our R mg with the results oflMaller et al.l (120061 ) at mass ratio m/M > 0.5 (Figure[ 
and the results of lStewart et all (120091 ) at mass ratio m/M > 0.3 and m/M > 0.6 (Figure[ 
We can see that both of Figures [25] and [26] show that our results are consistent with these 
previous results. 



Wetzel fc Tinker! (120121 ) found that from redshift z = 2 to now, it is around 27 percent 



of galaxies similar to our Milky Way that experienced a merger with mass ratio m/M > 0.1, 
and around 11 percent that experienced a merger with mass ratio m/M > 0.33. In our 
research, the galaxies in blue7 or blue8 subsample are similar to our Milky Way. If we take 
our results under assumption of t merge = txw,i at face value and make the strong assumption 
that the growth happens at a non-evolving rate, from redshift z = 2 to now (a period of 
~ 10 Gyr), we expect the galaxies in blue7 or blue8 subsample to merge by ~ 21 percent 
with mass rati o m/M > 0.1, and ~ 10 percent with mass ratio m/M > 0.33, which is close 



to the result of IWetzel fc Tinker! ( 12012! ) 
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5. Discussion 

We find that under the assumption of t merge = txw,i, the total fractional accretion 
rates onto red main galaxies are from [1.3 ± 0.7] to [3.7 ± 1.9] hm percent per Gyr depend- 
ing on stellar mass, and those onto blue main galaxies are from [0.6 ± 0.3] to [1.1 ± 0.6] 
/170 percent per Gyr. We find that at fixed stellar mass, the total fractional accretion rates 
onto red galaxies is greater than that onto blue galaxies. The total fractional accretion rate 
is a stronger function of primary mass for red galaxies than that for blue galaxies. We also 
find that more than 60 percent of the total fractional accretion rates are from major mergers 
with mass ratio between 1 : 1 and 1 : 3, and less than 15 percent of the total fractional 
growth rates are from minor mergers with mass ratio between 1 : 10 and 1 : 100. 

The first limitation of the imaging sample arises from the lack of spectroscopic informa- 
tion on the galaxies from the imaging sample. However, for galaxies from imaging sample 
with stellar masses Mi > 6 x 10 9 M , we estimate the real-space number densities rij 
from equation f l3"5"|) . So we measured rii, pi and w si {r p ) separately. But at stellar masses 
Mi < 6 x 10 9 /if 2 M , it is impossible to precisely measure the real-space number densities 
Hi, so that we cannot disentangle the clustering power from the number density for these 
small-mass galaxies from the imaging sample, and we only measure the products riiW S i{r p ) 
and PiW si (r p ) but not either rii, pi or w si {r p ) separately. 

The second limitation is removing galaxies from imaging sample with apparent magni- 
tudes r > 21.5 mag because of the limitation of lack of imaging information on the galaxies 
of SDSS (see Figure [5]). This cut of apparent magnitudes r > 21.5 mag will affect the minor 
mergers with mass ratio between 1 : 30 and 1 : 100 for the small stellar mass galaxies from 
spectroscopic samples redl~4 and bluel~4 at redshift z > 0.10, and the number densities 
of the galaxies in these spectroscopic samples decrease sharply at redshift z > 0.10 (see 
Figures |3] and @J, which will sharply reduce the effect of this cut. This cut will cause less 
than 5% error because the contribution of the minor mergers with mass ratio between 1 : 30 
and 1 : 100 is only < 5%. This assumption is good to take, because it will cause far below 
5% error. In order to not affect the minor mergers with mass ratio between 1 : 30 and 1 : 100 
at all, we need our imaging sample to be ~ 1 mag fainter than what we use now. 

We can see 'valleys' at the third to fifth data points (37.6 < r p < 94.4/ifQ 1 kpc) of the 
two minor-merger curves (10~ 1,5 < Mi/M s < 10 _1 and 10~ 2 < Mi/M s < 10~ L5 ) in each 
plot of Figures [TD] through [TSJ which seems like a kind of issue due to photometric biases 
or bad deblending. However we do not think so, because the photometric biases and bad 
deblending are very small at separations large than 15 arcsec (see Table [1]) and at 'valleys' 
(the third to fifth data points 37.6 < r p < QAAh^ kpc), the photometry correction is nearly 
zero (below 5%). Also, we do not think that the 'valleys' will make our conclusions invalid 
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even if that is an issue in minor merger, because from the above paragraph we know that 
only less than 15 percent of the total fractional accretion rates are from minor mergers with 
mass ratio between 1 : 10 and 1 : 100. 

If we assume that the growth happens at a non-evolving rate from redshift z = 1 to now 
(a period of « 8 h^Q Gyr), we expect the red galaxies to grow by about [10 ± 5] percent to 



merge 



t 



KW.i 



and the red 



[28 ± 14] percent depending on stellar mass under assumption of t 
L* galaxies (around red7 and red8) grow by about [20 ± 10] percent. The merger rate ma y 
have been different in the past, of course, higher or lower ( )Lin et al.ll2008l ; IChou et al.ll201ll ). 
If we assume that the growth of massive red galaxies (L* galaxies and above) is mainly from 
galaxy mergers, and also assume that the evolutio n of the galaxy merger rate per galaxy 
is proportional to (1 + z) +s 1-1 (ILotz et al.ll201ll ) from redshift z — 1 to now, we expect 
that the stellar mass density of the red massive galaxies (L* galaxies and above) increased 
about ~ 75 percent under the assumption of t merge = txw,i, or about ~ 40 percent under the 
assumption of t merge = tsT,i (see Table [3]), which are consistent with recent studies on the 
high-redshift Universe which find that the red sequence appears to grow in stellar mass over 



time by a factor of 50 percent to 1 percent from redshift z = 1 to now (|Bell et al. 



Willmer et al. 


2006: 


Blanton 


2006; 


Faber et al. 


2007 


2008; 


van Dokkum et al. 


2008 


; Saracco et al. 


20091). 



2007; Conselice et al. 200 



7; 



2004 



Cimatti et al. 



According to Sec tion HI our results are consistent with the previous clustering results 
of lZehavi et al] ( I2010l ). Our esti mated merger rates are consistent with the merger rates esti- 
mated by counting of close pairs ( Masiedi et al. I2OO8I Patton &; Atfieldll2008 ; Kitzbichler &; White 
20081 : iBundv et allbood : be Propris et al.lboiol ; iRobaina et aJboid ) . Our estimated merger 
rates are also consiste nt with the merger rates predicted in theories of galaxy format ion in a 
cosmological context jMaller et al.lbood : Istewart et alJliooalWetzel k Tinkerl lioia ) . 



However, we found that not all merger studies find such low values when we com- 
pared our results with the studies at higher redshift. The morphol o gical derivations of the 



merger fraction (e.g. , iDe Propris et al] 120071 : IConselice et al] 12007 



Lotz et al. 2006, 2008 



Conselice et al] 120091 ) tend to find values of f pa i r ~ 0.1 at 0.4 < z < 1.4, about 2 times 



higher than the results of our pair analysis (see Figure [24]). The discrepancy can be resolved 
easily if either (1) morphological signatures of merging last for many dynamical times (e.g., 
tidal tails) or (2) the very minor mergers (m/M < 0.25) inflate the merger rates (please 
note that our pair analysis in Figure [241 are only estimated from major mergers at mass ratio 
m/M > 0.25) or (3) morphological tools for finding mergers maybe find some systems which 
are not involved in mergers or (4) merger rates per galaxy at high redshift may be larger 
than those at low redshift. 



The total fractional accretion rates shown in the solid lines in Figure IT9l are upper limits 
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on the true fractional mass growth. There are two reasons: (1 ) we assume t merge ^ = tjxw,i 
for every pair. However it is not true for close pairs in high velocity dispersion, (2) we 
assume that the stellar mass growth of the central galaxies from the spectroscopic sample is 
equal to the stel lar mass of the galaxy from the imaging sample, see equation (|30p . However 
Lin et all (120041 ) found that up to 50 percent of the stars in the galaxies from the imaging 
sample could be stripped off before the merger with LRGs is complete. So our mass growth 
rate under the assumption of t merge = t KW)i is an upper limit on the growth by merging. 

We find that the accretion onto red and massive blue galaxies is dominated by mergers 
with red companions, and that onto small-mass blue galaxies, red and blue companions make 
comparable contributions, this is shown by Table |2] and Figures [T7] and [TBI So. most of the 



mass brought into red galaxies by mer ging is brought by "dry mergers" (jBell et al.ll2006b 
van Dokkumll2005l : iMasiedi et allboosl ). 



We find that all the contributions to growth decrease with decreasing stellar mass at the 
small-mass end for all of 18 spectroscopic subsamples. The contribution to growth decreases 
with decreasing Mi/M s since Mi/M s < 0.4. For all 18 subsamples, the curves essentially 



decrease to zero by Mi/M s < 0.01 for t 



merge 



t KW)i and by Mi/M s < 0.1 for t 



merge 



t 



BTa 



so calculation of the total amount of mass brought in by merger activities does not require 
consideration of galaxies from the imaging sample with Mi/M s < 0.01. 

From Figures [TUlthrough [TBI and Figure I2DI we find a sharp break at r p « 2.5 hjQ Mpc and 
a less-sharp transitio n at r p « 0.43 h? Mpc. These two transitions are also found in LRGs by 



Masjedi et all (120081 ), which can be exp 



of ga l axy clustering ([Peacock fc Smith 



ained in the context of the "ha 



2000 



Scoccimarro et al. 



2001 



o occupation" picture 



Berlind fc Weinberg 



20021 ; ICoorav & Sheth 



2002 



Zheng et al 



2005 



Watson et al 



2010, 



20121 ) : (1 ) the mergers 



at r p < 0.43 Mpc are the one-halo mergers (both of the two merging galaxies are inside 
one halo); (2) the mergers at r p > 2.5 Mpc are the two-halo mergers (the two merging 
galaxies are separately inside two nearby halos); (3) the mergers at 0.43 <r p < 2.5 Mpc 
are the mixed-halo mergers (some of the mergers are the one-halo mergers, the others are 
the two-halo mergers). So at r p « 0.43 hjQ Mpc the mergers transfer from the one-halo 
mergers to the mixed-halo mergers, and at r p « 2.5 hjQ Mpc the mergers transfer from the 
mixed-halo mergers to the two-halo mergers. These two transitions are clearer in Figure [TU] 
of mergers between two red galaxies than that in the other three figures. 
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Table 1. Spectroscopic Subsamples 
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Note. — Information for the galaxies in the spectroscopic subsamples. 

a Mean stellar mass of galaxies in the corresponding spectroscopic subsample. 

b Mean redshift of galaxies in the corresponding spectroscopic subsample. 

c Lower redshift limit used. 

d Upper redshift limit used. 

e We show adopted velocity dispersion a v for the red galaxies and adopted circular velocity 
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V c for the blue galaxies. 

f We show the angular separations of r p = 37.6/ifo kpc and r p = 59.6/ifQ 1 kpc which are 
the separations of the third and fouth data points in Figures [10] through [131 
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Table 2 
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Note. - 


Fractional mass g 


;rowth of Main Galaxy by merging Gyr 


, split by spectro- 



scopic subsample. 

a Measurements under the assumption t me rge,i = tsT,i- 

b Percent contribution of the blue galaxies from the imaging sample under the assumption 

Emerge, i ^BT,i- 

c Measurements under the assumption t merge ^ = txw,i- 

d Percent contribution of the blue galaxies from the imaging sample under the assumption 
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tmerge,i t>KW,i- 

e Errors are estimated including the error from merger time. 
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Table 3. Fractional mass growth measurements for massive red galaxies 

Subsample d(ln M) / d t | B7 EI d(ln M) / d * \ KV ^ 
[10- 3 h 70 Gyr- 1 ] [lO-^yoGyr" 1 ] 



red7 
red8 
red9 



31.8 ± 15. & 
33.8 ± 16.9 
40.5 ±20.3 



52.6 ±26.3 
58.1 ±29.1 
83.1 ±41.6 



> L* 
red7~9 



37.0 ± 18.5 70.5 ±35.3 



^ L* 
red8~9 



38.3 ±19.1 74.9 ±37.4 



LRG 

red9 



40.5 ±20.3 



83.1 ±41.6 



Note. — Fractional mass growth measurements for 
massive red galaxies (L* galaxies and above) under 
assumption of merger rates oc (1 + z )+ 3 0±11 ; S pHt by 
spectroscopic subsample. 

a Measurements under the assumption t merge .i = 

tBT,i- 

b Measurements under the assumption t merge ^ = 

c Errors are estimated including the error from 
merger time. 
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Spectroscopic Sample (3xl0 4 of 5xl0 5 ) 
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Fig. 1. — Sky distribution of the spectroscopic and imaging samples. For clarity, only a 
randomly chosen subsample of 3 x 10 4 points is shown in each case. 
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Fig. 2. — Distribution of the spectroscopic subsample in color and fT-correct-estimated 
stellar mass. The nearly horizontal line separates galaxies into red and blue; the vertical 
lines separate galaxies into 10 subsamples with different stellar mass. For clarity, only a 
randomly chosen subsample of 3 x 10 4 points is shown. 
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Fig. 3. — The number densities of each of the 9 red spectroscopic subsamples as a function of 
redshift. The vertical lines in each graph show the redshift limits used. The number density 
in red9 appears to rise with redshift, because we have removed galaxies with r < 14 mag (see 
text). 



-36 - 




0.05 0.10 0.15 0.05 0.10 0.15 0.05 

redshift 



0.10 0.15 



Fig. 4. — Same as Figure [3J but for the 9 blue spectroscopic subsamples. The number 
density in blue8 and blue9 appears to rise with redshift, because we have removed galaxies 
with r < 14 mag (see text). 
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Fig. 5. — Distribution of apparent r magnitude for the imaging sample. We use only galaxies 
from the imaging sample with 14 < r < 21.5 mag. 
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Fig. 6. — Distribution of [g — r] color for the imaging sample. We use only galaxies from the 
imaging sample with —0.5 < [g — r] < 2.0 mag. 




Fig. 7. — The difference between the i^-correct-estimated stellar mass and the stellar mass 
estimated by the Grid Method. The contours show the 68, 95 and 99-percent intervals. The 
greyscale reflects the number of galaxies in each bin. 
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Fig. 8. — The difference between the .fT-correct-estimated [° g — r] color and the [ ai g — 
01 r] color estimated by the Grid Method. The greyscale and the contours are similar to 
those in Figure [7J 
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Red galaxies from the Spectroscopic Sample 
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Fig. 9. — The relationship between a v , the stellar velocity dispersion and M, the stellar 
mass for red galaxies in the spectroscopic sample. The greyscale and the contours reflect the 
number of galaxies in each bin. The solid line is the linear fit, equation fl6]). 
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0.01 0.10 1.00 0.01 0.10 1.00 0.01 0.10 1.00 

r p (h^ Mpc) 

Fig. 10. — Projected two-dimensional cross-correlation functions PiW p (r p ) between all of 
the 9 red spectroscopic subsamples and all red imaging subsamples, scaled by r p for better 
illustration. The right vertical thin lines are corresponding to 55 arcsec at the median redshift 
of the spectroscopic galaxies, and the left vertical thick lines are corresponding to the median 
r 90 for the spectroscopic galaxies. The error-bars are from the jackknife error covariance 
matrix only. All the lines show the results after photometry correction: the thick solid line 
shows imaging galaxies with 1CT - 5 < Mi/M s < 10° (black diamonds are the result before 
photometry correction), the thick dashed line shows 1CT 1 < Mi/M s < 10~ a5 (black triangles 
are the result before photometry correction), the thin solid line shows 10~ L5 < Mi/M s < 10 _1 
(black squares are the result before photometry correction), and the thin dashed line shows 
< Mi/M s < 10~ L5 (black crosses are the result before photometry correction). The 
four curves are the fit lines(see text). The results before photometry correction are offset by 
12% of our interval to the left. 
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Fig. 11. — Same as Figure [101 but for the 9 red spectroscopic subsamples cross-correlated 
with all blue imaging subsamples. 
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Fig. 12. — Same as Figure [TOj but for the 9 blue spectroscopic subsamples cross-correlated 
with all red imaging subsamples. 
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Fig. 13. — Same as Figure [TOj but for the 9 blue spectroscopic subsamples cross-correlated 
with all blue imaging subsamples. 
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Fig. 14. — Recovered Petrosian flux to input Sersic flux as a function of the separation of the 
two galaxies in the pair. We show both LRG and small galaxy whose radius is only 80% of 
that of LRG. The vertical line shows the smallest separation in our research at z — .15 and 
r p = lA.9h^Q kpc. For LRG, on average there is an excess in the recovered flux of galaxies 
separated by less than 20 arcsec. 
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15. — Recovered Petrosian flux similar as Figure 1201 
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Fig. 16. — We show the major merger of LRG (IMasjedi et al.l 120061 ) . major merger of main 
galaxy and minor merger of main galaxy. Please note that this is only a sketch, the radius 
and separation may be much different. 
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Fig. 17. — The mean fractional accretion rate for the 9 red spectroscopic subsamples 
hfo Gyr -1 per dex. The thick lines are for mergers with red galaxies from the imaging 
sample and the thin lines are for mergers with blue galaxies from the imaging sample. The 
solid lines are the merger rate under assumption of t merge>i = tK\v,i, an d the dashed lines are 
for the merger rate under assumption of t merge ^ = tBT,i- 
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Same as Figure H7J but for the 9 blue spectroscopic subsamples. 
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Fig. 19. — The total fractional accretion rates for each of the 18 spectroscopic subsamples 
hro Gyr^ 1 integrating over all the galaxies from the imaging sample. The thick lines are for 
the red spectroscopic subsamples and the thin lines are for the blue spectroscopic subsamples. 
The solid lines are for the merger rate under assumption of t merge ^ = txwj, an d the dashed 
lines are for the merger rate under assumption of t merge ^ = tsT,i- 
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Fig. 20. — Projected correlation function w p (r p ) for the spectroscopic subsamples corre- 
sponding to mass-threshold samples as labeled, calculated as described in the text on small 



scales , combined with projected correlation function on intermediate scales from lZehavi et al. 



(201 0[). Plea s e not e that in order to compare these results easily, we offset the points of 



Zehavi et al.l ( 120101 ) by 12% of our interval to the right. 
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Fig. 21. — Projected correlation function w p (r p ) similar as Figure EUJ There is no offset in 
this graph. The triangle points are our w, 
points ar e the w p (r p ) a s sumin g; pj = 1 and fj 



p\T p ) with correction of fiber collisions, the diamond 
= 1, the solid lines with thin error bars are the 



w p (r p ) of IZehavi et al.l ( 120101 ) and the dashed lines are the extension lines of IZehavi et al. 
(120101 ) described in the text. Please note that in order to compare these results easily, we 
only display our first ten data points. 
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Fig. 22. — Projected correlation function w p (r p ) similar as Figure |2T| but only display the last 
three plots. The black crosses in the second and third plots are our result before photometry 
correction which are offseted by 12% of our interval to the left. Please note that the three 
extension da s hed li nes are fitting from the first five, first six and first seven data points of 
Zehavi et all feoioh . 
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Fig. 23. — The comparison of our major merger rate MF at 0.03 < z < 0.15 under assumption 
of t merge = txw,i with the major merger rate MF of IBundy et al.l (120091 ) at 0.4 < z < 0.7. 
Both of our mass ratios are m/M > 0.25. The open diamonds on t he thick lines are ou r 
results, and the open triangles on the thin lines are the results of IBundy et al.l ( 120091 ). 
The dashed lines show the results of excluding the approximate fraction of dry E/S0-E/S0 
mergers, and the solid lines indicate the observed merger rate for all galaxies determined. 
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Fig. 24. — The comparis on of our f 



7.1-28.6 kpc 



pair 



at 0.03 < z < 0.15 with the 



„7.1-28.6/if 1 kpc 
/ pair 



of th ree recent results (IPatton k Atfieldl 120081 : Kitzbichler k White! 120081 : iBundv et al. 
2009h at 0.4 < z < 0.7. The open diamonds on the thick solid lines are the results of 



Bundy et al.l (120091 ). and the open triangles o n the thick dashed lines are our results. The 



thin horizontal solid lines show the res ults of IPatton k Atfieldl (|2008T) . and the thin hori- 
zontal dashed lines show the results of Kitzbichler k White! ( 20081 ). Please note that the 
Z p 7 al~s 28 ' 6h7 ~° lkpC of IPatton k Atfieldl J2008h and iKitzbichler k White! J2008h are calculated by 

fpair 
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Fig. 25. — The comparison of our merger rates (per gala xy hm Gyr 1 ) unde r assumption of 



t 



merge 



tK\v,i a t 0.03 < z < 0.15 with the merger rates of lMaller et al.l (120061 ) at < z < 0.5. 



The open triangles on the dashed lines are our results, and the thin dotted lines are the best- 



fit result calculated from equation (5) of iMa" 



ler et al. ( 2006 ). Both of our mass ratios are 



m/M > 0.5. Please note that the results of 
Q m = 0.4 cosmology with a 8 

h = 0.65, a baryon content fib = 0.047, and a spectral index n 



Mailer et al.l (120061 ) are estimated from a flat 
0.8, a Hubble constant H = 100 /ikms" 1 Mpc -1 with 

0.93. This difference of the 



two assumptions might lead the difference of the two slopes. 
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Fig. 26. — The comparison of our merger rates (per galaxy Gyr -1 ) under assumption 
of t merge = txw,i at 0.03 < z < 0.15 with the estim ated merger rates calculated by the 
'Merger Rate Fitting Function' of iStewart et al.l (120091 ) at z — 0.1. The open diamonds on 
the thick lines a re our results, and th e horizontal thin lines are the best-fit results calculated 
from Table 1 of IStewart et al.l (120091 ). The dashed lines show the results with mass ratios 
m/M > 0.3, and the solid lines show the results with mass ratios m/M > 0.6. Please 
note that the results of IStewart et al.l (120091 ) are estimated from a flat Q m = 1 — Q\ = 0.3 
cosmology with a Hubble constant Hq = 100 /ikms -1 Mpc^ 1 and h = 0.7. 



